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Foreword 


The  first  three  years  (04/23/2008  -  06/22/2011)  of  this  program  involved  the  joint  efforts 
of  two  groups,  the  synthesis  group  with  Prof.  Karl  Christe  as  PI  and  the  simulation  group 
with  Prof.  Priya  Vashista  as  the  PI.  On  June  22,  2011,  DTRA  granted  a  two  year 
extension  to  the  synthesis  effort  by  the  Christe  group,  and  awarded  Dr.  Ken-ichi  Nomura 
a  2-year  Young  Investigator  Award  to  continue  the  exploration  of  the  computational  part 
of  the  project.  This  report  summarizes  the  results  obtained  by  both  groups  during  the 
original  program  and  the  results  obtained  by  the  Christe  group  under  the  synthesis 
extension  program. 


Objectives 

Chemical  and  biological  agents  (CBAs)  represent  a  serious  threat  to  our  homeland 
security.  Once  these  agents  become  airborne  over  large  areas,  they  become  very  difficult 
to  destroy.  Consequently,  it  is  desirable  to  either  destroy  these  agents  before  they  can 
escape  from  their  storage  or  production  facilities  or,  even  better,  to  destroy  them 
simultaneously  with  the  entire  facility.  One  can  envision  a  combined  weapon  system  in 
which  tailored  explosives  would  destroy  the  storage  or  production  facility,  while 
secondary  chemical  reagents  would  neutralize  the  CBAs.  Therefore,  the  objectives  of  this 
proposal  were  the  development  of  powerful  novel  chemical  reagents  for  the  destruction 
of  biologically  active  materials  and  a  simulation  of  their  reactions  on  a  multimillion  atom 
scale  with  quantum  mechanical  accuracy  to  study  the  stability  of  the  reactive  species 
under  high-temperature  and  high-pressure  fireball  conditions. 

More  specifically,  we  have  been  synthesizing  and  characterizing  novel  halogen  oxidizers 
for  formulation  with  high-performance  explosives  (HMX  or  RDX)  to  produce  over¬ 
oxidized  explosives  for  destruction  of  chemical  &  biological  agents.  Multimillion-atom 
molecular  dynamics  simulations  with  quantum  mechanical  accuracy  were  used  to  study 
the  stability  of  the  halogen  species  in  high-temperature  and  high-pressure  fireballs.  Ideal 
candidates  for  these  over-oxidized  explosives  were  mixtures  of  HMX  or  RDX  with  halogen 
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oxides  in  their  highest  oxidation  states.  Of  these  halogen  oxides,  the  unknown  10^  originally 
was  our  preferred  choice.  It  was  predicted  to  be  thermally  stable  and  insensitive  and  to 

3 

possess  a  density  in  excess  of  4.5  g/cm  .  Therefore,  the  premier  goal  of  our  proposed  study 

was  the  synthesis  and  characterization  of  this  compound.  However,  other  sources  for  positive 

+  + 

halogen  compounds,  such  as  molecular  I206,  lO  F,  R^O,  and  IF  O^,  IF^O,  ionic  IF  O  ,  IO^  , 

-  -  -  2-2-  +  - 

IF  O  ,  IF^O  ,  IF  O^ ,  IF502  ,  IF_0  and  IF^  IO^  were  also  good  target  compounds. 

The  original  program  plan  included  the  following  milestones. 

Year  1:  Synthesis  and  characterization  of  I2O7  precursors,  theoretical  prediction  of 
properties  of  I2O7  and  modeling  of  iodine  oxide/HMX  using  the  known  compounds  I2O5 
and  I206. 

Year  2:  Synthesis  and  characterization  of  I2O6  and  modeling  of  the  I2O7  and  I2O6 
systems. 

Year  3:  Development  of  Scale-up  synthesis  and  formulation  testing  of  I2O6  and  modeling 
of  the  complete  system  including  CBAs. 

The  goals  of  the  2-year  Synthesis  Extension  were:  Preparation  of  larger  amounts  of  I2O6 
and  their  distribution  to  other  laboratories  for  advanced  formulation  studies  and  testing; 
crystal  structure  determination  of  I2O7;  and  the  development  of  coating  techniques  for 
I2O6  to  improve  its  hydrolytic  stability. 

Summary  of  Efforts 

Synthesis  Task 

Excellent  progress  was  made  during  this  program.  We  have  successfully  synthesized  the 
previously  unknown  I2O7  and  have  characterized  it.  The  only  missing  piece  of 
information  is  a  crystal  structure.  We  have  shown  that  it  is  of  insufficient  thermal 
stability,  but  that  the  known  I2O6  is  an  excellent  substitute.  We  have  drastically  improved 
the  literature  synthesis  for  I2O6,  have  scaled  up  its  production  to  the  100  g  level,  have 
prepared  several  hundred  gram  of  I2O6,  and  have  provided  Prof.  Yashistha  (USC),  Prof. 
Dreizin  (NJ  Inst,  of  Technology),  Prof.  Zachariah  (U  of  Maryland),  Dr.  Higa  (NAY AIR, 
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China  Lake),  Dr.  Lightstone  (Indian  Head),  Dr.  Liu  (NIST)  and  Dr.  Kolesnikov 
(Oakridge)  with  samples  for  further  evaluation.  Explosives  testing  at  NAVAIR  and 
Indian  Head,  powder  diffraction  measurements  at  NIST  and  neutron  diffraction  and 
inelastic  scattering  experiments  at  Oakridge  have  been  highly  successful.  At  the 
suggestion  of  Dr.  Su  Peiris,  we  have  also  carried  out  coating  studies  for  I2O6  which 
resulted  in  an  improvement  of  its  hydrolytic  stability. 

Simulation  Task 

We  have  performed  first-principle  quantum-mechanical  calculations  of  I2O5  and  I2O6 
crystals  to  investigate  their  structural  and  lattice-dynamic  properties.  We  have  also 
calculated  their  equations  of  states  and  thermodynamic  properties  at  high  pressures  and 
temperatures.  The  results  were  summarized  in  a  manuscript  entitled:  “Phonon  Density  of 
States  of  Iodine  Oxides  -  I2O5  and  I2O6:  A  Joint  Study  Using  Neutron  Scattering  and 
First-principles  Calculations,”  which  will  be  submitted  to  Chem.  Phys.  Letters  for 
publication. 


Accomplishments/New  Findings 

Synthesis  Task 

The  synthesis  of  the  previously  unknown  I2O7  was  achieved  by  the  low-temperature 
dehydration  of  H5IO6  with  SO3  (eq.l). 

H  SO 

2  H5I06  +  5  SO3  — 2  4  ►  I207  +  5  H2SO4  (1) 

The  I2O7  was  identified  by  the  observed  material  balance  and  its  Raman  spectrum  which 
differs  significantly  from  that  of  I2O6  (see  Fig.  1).  The  observed  Raman  spectrum  of  I2O7 
was  also  compared  to  that  for  the  free  I2O7  molecule  calculated  for  us  at  the 
B3LYP/DZYP  level  by  Prof.  David  Dixon  at  the  University  of  Alabama.  The  comparison 
shows  that  the  observed  spectrum  is  much  more  complex  than  that  calculated  for  the  free 
gaseous  molecule,  and  that  I2O7  must  possess  an  oxygen-bridged  polymeric  structure. 
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Figure  1.  Raman  spectra  of  I2O7  and  I2O6 

Therefore,  a  crystal  structure  is  needed  to  establish  the  exact  structure  of  solid  I2O7. 
Efforts  to  grow  single  crystals  of  I2O7  were  made  but  were  impeded  by  its  low  solubility 
in  all  the  tested  solvents.  The  I2O7  is  a  marginally  stable  yellow  solid  which  starts  to  lose 
dioxygen  at  room  temperature  and  decomposes  rapidly  at  60  °C  to  give  I2O6  (eq.  2). 

2  I207  - ►  2  I206  +  02  (2) 
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Because  of  its  low  thermal  stability,  I2O7  is  not  a  good  candidate  for  practical 
applications.  However,  our  Thermo  Gravimetric  Analysis  (TGA)  and  Differential 
Scanning  Calorimetry  (DSC)  studies  of  I2O6  showed  that  this  compound  is  thermally 
stable  to  about  187  °C,  where  it  decomposes  quantitatively  to  I2O5  and  dioxygen  (eq.  3). 

2  I206  - ►  2  I205  +  02  (3) 

At  about  400  °C,  the  I2O5  then  decomposes  to  iodine  and  oxygen.  A  special  differential 
scanning  calorimeter  was  built  (see  Fig.  2)  and  used  to  measure  (see  Fig.  3)  the  heat  of 
formation  of  I206.  The  DSC  curve  showed  an  exotherm  with  an  onset  of  179  °C  and  a 
peak  of  197  °C  for  the  decomposition  of  I206  to  I2O5  and  dioxygen,  and  an  endotherm 
with  an  onset  of  393  °C  and  a  peak  of  424  °C  for  the  decomposition  of  I2O5  to  the 
elements.  From  these  data,  the  heat  of  formation  of  I206  was  found  to  be  -90.2  kJ/mol. 


Figure  2.  Differential  Scanning  Calorimeter 
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Temperature  (*C) 

1)  Hold  for  10  mm  at  30.00 C  3)  Cool  from  700  00 to  30.00"C  at  10  OO’C/min 

2)  Heat  from  SO.OO'C  to  700.00^  at  10.00*C/min 


Figure  3.  DSC  Curve  of  I2O6 


This  heat  of  formation,  combined  with  its  density  of  5.40  g/cm3  (T.  Kraft  and  M.  Jansen, 
J.  Am.  Chem.  Soc.,  1995,  117,  6795)  and  its  ability  to  generate  spore  killing  positive 
halogen  compounds,  makes  I2O6  an  outstanding  candidate  for  the  DTRA  applications. 

The  crystal  structure  of  I2O6  is  well  known  (T.  Kraft  and  M.  Jansen,  J.  Am.  Chem.  Soc., 
1995, 117,  6795),  but  the  published  synthesis  from  H5IO6  and  cone.  H2SO4  had  been  very 
time  consuming.  We  have,  therefore,  worked  on  an  improved  synthesis  of  I2O6.  Instead 
of  going  through  the  decomposition  of  an  intermediate  I2O7,  we  have  used  equimolar 
amounts  of  iodic  and  periodic  acids  instead.  In  this  way,  we  have  shortened  the  reaction 
time  from  the  previous  one  month  to  one  day.  We  have  scaled  up  the  reaction  to  a  100  g 
level  and  have  prepared  several  hundred  g  of  I2O6.  We  have  provided  Prof.  Vashistha 
(USC)  with  large  samples  for  the  neutron  diffraction  studies.  We  have  also  sent  samples 
of  I2O6  to  Prof.  Dreizin  (NJ  Inst,  of  Technology),  Prof.  Zachariah  (U  of  Maryland),  Dr. 
Higa  (NAVAIR,  China  Lake),  Dr.  Lightstone  (Indian  Head),  Dr.  Liu  (NIST)  and  Dr. 
Kolesnikov  (Oakridge)  with  samples  for  further  evaluation.  We  have  also  prepared  a 
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sample  of  hydrogen-free  I2O5  by  pyrolysis  of  I2O6  for  a  neutron  diffraction  study  by  Dr. 
Kolesnikov  at  Oak  Ridge. 

Dr.  Higa  at  NAVAIR  has  obtained  very  interesting  results.  While  neat  I2O6  reacted  only 
slowly  with  80  nm  A1  powder,  the  same  mixture  with  an  AI/M0O3  bootstrap  gave 
deflection  rates  about  twice  as  large  as  those  of  an  Al(320  nm)/Bi203  standard. 
Furthermore,  the  AI/I2O6  composite  was  stable  on  storage  for  more  than  10  months.  Very 
importantly,  in  all  of  these  tests,  the  formation  of  elemental  iodine  was  visually  observed. 
Similar  tests  with  I2O5  in  place  of  I2O6  resulted  only  in  half  of  the  deflection  rates  and 
demonstrate  that  I2O6  is  clearly  superior  to  I2O5.  A  patent  disclosure  on  these  I2O6  based 
explosives  has  been  submitted  by  Dr.  Higa  through  NAVAIR. 

Spore  killing  tests  with  Al/I206  formulations  were  carried  out  at  NSWC  Indian  Head  by 
Jillian  Home  and  Dr.  Jim  Lightstone.  They  found  that  the  I2O6  formulations  resulted  in 
greatly  improved  temperature/time  plots  both  in  the  detonation  chamber  and  the  exhaust 
chamber.  Due  to  the  unexpectedly  high  pressures  obtained  in  these  tests,  the  filters  were 
blown  out  and  they  did  not  obtain  data  on  the  spore  killing  efficiency  of  these 
formulations.  The  tests  must  be  repeated. 

Since  I2O6  is  somewhat  hygroscopic,  Dr.  Su  Peiris  suggested  carrying  out  coating  studies 
to  improve  its  hydrolytic  stability.  In  preliminary  experiments,  it  was  attempted  to  coat 
I2O6  with  perfluorodecanoic,  perfluorododecanoic  and  perfluorooctadecanoic 
(perfluorostearic)  acid.  For  the  experiments,  samples  of  about  1  g  of  finely  powdered 
I2O6  was  stirred  in  about  10  mL  of  10  weight  %  solutions  of  the  perfluorinated  acid  in 
trifluoroacetic  acid.  The  trifluoroacetic  acid  was  then  evaporated  in  a  vacuum  so  that  the 
dissolved  perfluorinated  acid  resulted  in  a  coating  of  the  I2O6  particles.  The  obtained 
yellow  powders  were  inspected  by  Raman  spectroscopy.  In  all  three  cases,  the 
characteristic  bands  of  I2O6  together  with  the  much  weaker  bands  of  the  perfluorinated 
acid  were  present  in  the  spectrum.  The  I2O6  samples  were  then  exposed  to  the  air  and 
their  decomposition  followed  by  Raman  spectroscopy.  In  all  three  cases  the  samples  were 
less  air  sensitive  than  a  sample  of  untreated  I2O6.  However,  after  an  air  exposure  of  about 
4  hours,  Raman  bands  characteristic  for  I2O5,  HIO3  and  H5IO6  started  to  appear.  This  is 
clearly  an  indication  that  it  is  possible  to  make  I2O6  more  resistant  to  atmospheric 
moisture  by  coating  it  with  perfluorinated  compounds;  however,  more  experimental  work 
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in  this  area  is  needed.  Future  experiments  would  involve  the  use  of  perfluorinated  or 
partially  fluorinated  alkohols  such  as  ultrahigh  molecular  weight  2- 
(perfluoroalkyl)ethanol,  perfluorododecan-l-ol,  perfluorotetradecan-l-ol  and 
perfluorohexadecan-l-ol  as  well  as  solutions  of  amorphous  fluoropolymer  (Teflon  AF, 
“liquid  Teflon”). 

Simulation  Task 

I.  Structural  Properties 

The  I2O5  crystal  has  the  monoclinic  structure  and  space  group  P2  ]/c  with  four  molecules 
in  the  unit  cell  (see  Fig.  4),  whereas  the  I2O6  crystal  has  the  triclinic  structure  and  belongs 
to  space  group  P\  (see  Fig.  5). 


Figure  4.  Crystal  structure  of  I2O5. 


C 


Figure  5.  Crystal  structure  of  I2O6. 
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The  volumes  calculated  by  DFT  without  van  der  Waals  (VDW)  correction  are  far  larger 
than  the  experimental  equilibrium  volumes  for  the  two  iodine  oxide  crystals.  Namely,  the 
volume  difference  is  9.9%  and  11.7%  for  I2O5  and  I2O6,  respectively.  Such  discrepancy  is 
common  in  molecular  crystals  and  has  been  ascribed  to  the  inadequacy  of  exchange 
correlation  functional  in  dealing  with  VDW  interaction.  After  including  the  VDW 
correction,  the  calculated  volume  agrees  well  with  the  experiment  data. 

Calculation  without  VDW  correction  usually  overestimates  the  bond  lengths  about  2-3%. 
In  some  bonds,  the  overestimation  can  be  as  large  as  5%.  Although  VDW  correction 
definitely  decreases  the  volume,  its  effect  on  the  bond  length  is  not  trivial.  In  I2O5,  VDW 
correction  in  fact  increases  the  bond  lengths  slightly  except  for  one  bond.  The  change  of 
bond  length  caused  by  VDW  correction  is  usually  smaller  than  1%.  For  I2O6,  the  effect  of 
VDW  on  bond  lengths  depends  on  specific  iodine  atoms.  The  bonds  from  one  I  atom 
increase  lengths,  whereas  bonds  from  another  I  atom  decrease  lengths  by  VDW 
correction.  However,  the  change  of  the  bond  length  is  usually  small. 

For  intermolecular  distances,  VDW  correction  always  decreases  them.  In  general,  VDW 
correction  makes  these  interatomic  distances  more  consistent  with  the  experimental 
value. 

The  effect  of  pressure  on  bond  lengths  is  found  to  be  similar  to  that  of  VDW  correction. 
If  the  bond  length  increases/decreases  with  VDW  correction,  it  increases/decreases  with 
pressure.  Specifically,  all  the  bond  length  of  I2O5  but  one  increases  with  the  pressure.  For 
I2O6,  the  bonds  from  one  I  atom  increase  the  length  and  bonds  from  another  I  atom 
decrease  the  length  with  the  pressure.  Therefore,  to  some  extent,  the  effect  of  the  VDW 
correction  can  be  viewed  as  the  pressure  effect. 

II.  Vibrational  Properties 

The  phonon  density  of  states  (DOS)  of  I2O5  (Fig.  6)  has  a  small  gap  located  between  425 
and  439  cm'1.  The  detailed  analysis  of  the  vibrational  modes  (Fig.  7)  indicates  that  the 
stretching  modes  and  bending  modes  are  separated  by  this  gap.  Above  this  gap,  all  the 
vibrational  modes  are  stretching  modes.  There  is  another  gap  at  a  frequency  of  about  700 
cm'1,  where  the  phonon  DOS  of  I2O5  is  considerably  smaller  than  the  rest.  Above  this 
frequency  are  the  terminal  1=0  stretching  modes  and  below  it  are  the  bridging  I-O-I 
stretching  modes.  The  frequencies  of  the  stretching  modes  are  sensitive  to  the  bond 
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lengths.  The  calculation  overestimates  the  bond  lengths.  This  leads  to  an  overall 
underestimation  of  the  frequencies  of  the  stretching  modes.  In  contrast,  the  calculated 
frequencies  agree  well  with  the  experimentally  observed  ones  for  the  bending  and  the 
lattice  modes. 


Figure  6.  Partial  phonon  density  of  states  of  I2O5. 

The  Grlineisen  parameter,  which  is  defined  as  yt  =-d\ncoJd\nV ,  describes  the  volume 
dependence  of  the  frequency.  Since  the  bond  lengths  can  increase  slightly  with  the 
pressure,  some  of  the  stretching  modes  have  small  but  negative  mode  Griineisen 
parameters.  The  Griineisen  parameters  of  the  lattice  modes  of  I2O5  are  usually  larger  than 
1  and  also  significantly  larger  than  those  of  other  modes.  This  reflects  the  fact  that  inter- 
molecular  distances  are  greatly  reduced  by  pressure. 

Compared  to  solid  I2O5, 12O6  has  different  vibrational  properties.  For  example,  the  single 
bond  stretching  modes  are  not  separated  clearly  from  the  deformation  modes. 
Furthermore,  the  highest  two  stretching  modes  due  to  doubly  bonded  oxygen  have 
frequencies  far  larger  than  the  other  modes,  which  lead  to  a  broad  gap  in  the  phonon  DOS 
(see  Figs.  7  and  8).  Again,  these  modes  must  be  assigned  to  terminal  1-0  groups.  These 
are  also  the  only  two  modes  with  small  negative  Griineisen  parameters.  As  in  I2O5,  the 
Griineisen  parameters  of  the  lattice  modes  are  usually  larger  than  1,  which  is  also  the  case 
for  most  other  modes.  Since  the  calculation  also  overestimates  the  bond  lengths  of  I2O6, 
we  believe  that  the  frequencies  of  the  stretching  modes  are  also  underestimated  in  our 
calculation. 
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Figure  7.  Vibrational  modes  of  I2O5. 
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Figure  9.  Vibrational  modes  of  I2O6. 


III.  Equation  of  States  and  Thermodynamic  Properties 

The  equation  of  states  of  I2O5  and  I2O6  are  shown  in  Figs.  10  and  11.  The  calculation 
without  VDW  correction  overestimates  the  volume  significantly.  However  at  high 
pressures,  the  effect  of  VDW  on  volume  becomes  less  for  both  crystals.  The  volume 
differences  between  calculations  with  and  without  VDW  correction  for  I2O5  and  I2O6 
at  0  GPa  are  12.0%  and  13.2%,  respectively.  At  5  GPa,  these  differences  decrease  for 
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I2O5  and  I2O6  to  4.5%  and  5.7%,  respectively.  The  tendency  that  the  VDW  effect 
decreases  with  pressure  is  consistent  with  previous  results  for  several  other 
molecular  crystals. 


Figure  10.  Equation  of  states  of  I2O5. 


Figure  1 1 .  Equation  of  states  of  I2O6. 

Vibrational  contributions  increase  the  equilibrium  volume  noticeably.  For  I2O5  and  I2O6, 
the  effect  of  zero-point  motion  and  room  temperature  shifts  the  equilibrium  volume  by 
2.2%  and  0.9%,  respectively.  Interestingly,  these  values  are  very  similar  to  those  of 
materials  such  as  MgO  and  Mg2Si04,  although  the  bulk  moduli  of  MgO  and  Mg2Si04  (at 
about  170  GPa)  are  far  larger  than  those  of  iodine  oxide  crystals  (at  about  30  GPa).  The 
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calculated  equilibrium  volume  including  VDW  and  phonon  contributions  agrees  well 
with  the  experimental  data  (with  differences  smaller  than  0.5%).  The  resulting  equation 
of  states  provides  reliable  volumes  at  high  pressure  and  temperature. 

Thermodynamic  properties  of  two  crystals  including  thermal  expansion,  heat  capacity, 
adiabatic  bulk  modulus  and  thermal  Griineisen  parameters  are  shown  in  Figs.  12  and  13. 


Figure  12.  Thermodynamic  properties  of  I2O5:  (a)  thermal  expansion,  (b)  heat  capacity  at 
constant  pressure  Cp  and  volume  CV,  (c)  adiabatic  bulk  modulus,  and  (d)  Griineisen 
parameters. 
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Figure  13.  Thermodynamic  properties  of  I2O6:  (a)  thermal  expansion,  (b)  heat  capacity  at 
constant  pressure  Cp  and  volume  Cy,  (c)  adiabatic  bulk  modulus,  and  (d)  Griineisen 
parameters. 


The  thermal  expansion  coefficient  increases  rapidly  at  low  temperatures  up  to  about 
250K  and  then  increases  slightly  with  increasing  temperature.  The  bulk  modulus 
decreases  almost  linearly  with  temperature.  The  thermal  Griineisen  parameter,  which  is 
defined  as  yth  =  aKTV / Cv,  is  the  weighted  average  of  the  mode  Griineisen  parameters.  At 
low  temperatures,  the  thermal  Griineisen  parameter  is  determined  mostly  by  the  low 
frequency  modes.  With  increasing  temperature,  all  modes  tend  to  contribute  evenly  and, 
hence,  the  thermal  Griineisen  parameter  is  close  to  the  average  value  of  the  mode 
Griineisen  parameters.  Therefore  the  thermal  Griineisen  parameter  is  almost  temperature 
independent  at  sufficiently  high  temperatures.  Low  frequency  modes  (lattice  modes)  have 
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larger  Grtineisen  parameters  than  other  modes,  and  this  explains  why  the  thermal 
Griineisen  parameters  decrease  with  temperature. 

The  heat  capacity  at  constant  pressure,  Cp,  and  constant  volume,  CV,  are  related  by 
CP  =  (1  +  aythT)Cv .  Heat  capacity  increases  rapidly  below  about  250K  and  tends  to 
saturate  at  high  temperatures.  Since  we  use  a  quasi-harmonic  approximation  in  the 
calculation,  Cy  is  subjected  to  the  law  of  Dulong  and  Petit.  However,  Cp  can  be 
significantly  larger  than  the  Dulong  and  Petit  limit. 


IV.  Neutron  Scattering  Results 

Neutron  diffraction  and  inelastic  scattering  experiments  were  done  in  collaboration  with 
the  Spallation  Neutron  Source  at  Oak  Ridge  (Dr.  Alexander  Kolesnikov)  and  the  NIST 
Reactor  (Dr.  Yun  Liu).  Diffraction  measurements  at  low  temperature  were  completed. 
Analysis  of  the  inelastic  neutron  measurements  has  also  been  completed.  The  results 
from  these  studies  have  been  summarized  in  the  following  manuscript  for  publication. 
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Phonon  Density  of  States  of  Iodine  Oxides  -  I2O5  and  I2O6: 

A  Joint  Study  Using  Neutron  Scattering  and  First-principles  Calculations 

A.  I.  Kolesnikov,1  D.  L.  Abernathy,1  A.  Huq,1  J.  P.  Hodges,1  G.  E.  Granroth,1  C.-K. 
Loong,2  K.  O.  Christe,3  R.  Haiges,3  R.  K.  Kalia,4  A.  Nakano,4  K.  Nomura, 4  P. 

Vashishta,4  and  Z.-Q.  Wu4 

Neutron  Scattering  Sciences  Division,  Oak  Ridge  National  Laboratory,  Oak  Ridge,  TN 

37831-6473,  USA 

2 

School  of  Physics  and  Engineering,  Sun  Yat-Sen  University,  Guangzhou,  510275,  China 
Loker  Research  Institute  and  Department  of  Chemistry,  University  of  Southern 
California,  Los  Angeles,  CA  90089-1661,  USA 
4  Collaborator  for  Advanced  Computing  and  Simulations,  University  of  Southern 
California,  Los  Angeles,  CA  90089-0242,  USA 

Iodine  oxides  I2O 'y  (y  =  4,  5  and  6)  crystallize  into  atypical  structures  that  fall  in  between 
the  molecular-  and  framework-base  types  and  exhibit  high  reactivity  in  ambient 
environment.  We  performed  neutron,  Raman  and  IR  measurements  of  the  crystal 
structures  and  lattice  dynamics  of  newly  synthesized,  pure  I2O5  and  I2O6  samples.  The 
molecular  units  within  the  crystals,  their  connectivity  with  neighboring  units  to  form 
chain-like  or  framework-like  structures,  and  the  intra-and  inter-molecular  vibrations  and 
phonon  densities  of  states  were  characterized.  The  experimental  information  enables  the 
systematic  optimization  of  the  structures  and  dynamics  of  the  three  oxides,  I2O4, 12O5  and 
I2O6,  by  first-principles  calculations.  The  van  der  Waals  interaction  as  a  crucial  factor  in 
the  optimization  process  on  top  of  the  density  functional  theory  was  revealed  through  a 
quantitative  comparison  of  the  calculated  crystal  structures  and  phonon  densities  of  states 
with  experiments.  These  results  permit  the  calculations  of  the  equations  of  states,  bulk 
modulus,  thermal  Griineisen  parameters,  and  other  physical  quantities  that  are  useful  for 
further  understanding  of  the  thermodynamic  properties  of  iodine  oxides  under  extreme 
conditions. 
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Iodine  oxides  play  a  critical  role  in  many  important  physiochemical  phenomena  and, 
once  their  basic  thermodynamic  and  chemical  properties  are  understood,  may  lend 
themselves  to  a  series  of  overarching  technological  applications.  Firstly,  iodine  oxides, 
generated  by  photolysis  of  biogenic  iodocarbons  emitted  from  marine  algae,  interacts 
with  O3  under  ultraviolet  radiation  to  produce  aerosol  and  cloud  condensation  nuclei, 
hence  contributing  to  climate  change  [1,  2].  Secondly,  iodine  oxides  are  versatile  and 
effective  oxidants.  For  example,  I2O5  is  known  to  the  oxidation  of  alcohol,  the  nine- 
membered  amide  and  cycloalkan[Z?]indoles  whereby  useful  products  can  be  synthesized 
in  atom-efficient,  chemoselective,  and  environmental  friendly  way  [3,  4].  Finally,  I2O6 
and  I2O7  are  predicted  to  be  effective  neutralizing  reagents  that  defeat  the  functionality  of 
chemical  and  biological  agents. 

Notwithstanding  a  long  history  of  research,  experimental  characterization  of  iodine 
oxides  has  been  limited  mainly  to  I2O4,  I2O5,  and  a  few  intermediate  compounds 
associated  with  interactions  with  water  and  sulfuric  acid.  Ab  initio  calculations  have  been 
performed  in  various  molecular  iodine  oxides,  INO3,  HOI  species  [5,  6],  but  apparently 
no  rigorous  calculation  of  the  structural  and  thermodynamic  properties  of  bulk  I2O4, 12O5 
and  I2O6  has  appeared  in  the  literature.  Among  the  hOy  (y  =  4-7)  materials,  the  crystal 
structures  of  the  y  =  4,  5  and  6  members  are  found  to  conform  to  the  trend  of  changing 
from  polymeric-like  to  framework-like  structures  with  increasing  oxygen  content. 
Consequently,  the  thermodynamic  properties  of  I2O6,  particularly  under  high-pressure 
and  high-temperature  conditions,  are  expected  to  differ  from  those  of  I2O4  and  I2O5.  To 
this  end,  the  properties  of  I2O6  are  of  most  interest,  but  chemical  instability  of  this 
material  has  precluded  the  synthesis  of  bulk  I206,  let  alone  experimental  characterization. 
In  this  paper,  we  present  a  joint  neutron-diffraction/spectroscopic  plus  Raman/IR 
measurements  and  first-principles  calculation  of  the  structure,  phonon  density  of  states, 
and  thermodynamic  properties  of  bulk,  anhydrous  ^Oy.  The  theoretical  results  are 
compared  with  the  experimental  neutron,  Raman,  and  IR  data  in  the  literature  as  well  as 
from  present  measurements.  With  both  the  structure  and  lattice  dynamics  optimized,  we 
carried  out  calculations  of  the  equations  of  states  and  other  thermodynamic  properties  of 
the  three  iodine  oxides. 
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The  previously  in  the  literature  described  methods  for  the  synthesis  of  I2O6  have 
serious  drawbacks.  The  thermal  decomposition  of  H5IO6  in  a  vacuum  [7]  is  difficult  to 
control  and  produces  I2O6  that  is  contaminated  with  side  products.  Furthermore,  this 
preparation  is  only  applicable  for  very  small  amounts  of  material.  While  the  dehydration 
of  H5IO6  in  concentrated  H2SO4  at  70  °C  was  reported  to  proceed  only  very  slowly 
during  the  course  of  one  month  [8],  dehydration  of  a  mixture  of  H5IO6  and  HIO3  with 
65%  oleum  resulted  in  impurities  of  (102)28207  [9].  During  the  course  of  our  study,  it 
was  found  that  the  dehydration  of  a  mixture  of  H5IO6  and  HIO3  with  a  stoichiometric 
amount  of  SO3  in  100%  sulfuric  acid  at  elevated  temperature  yields  pure  I2O6  on  a 
several  gram  scale  within  a  few  hours: 

H5I06  +  HI03  +  3S03  H2S°4  >  I206  +  3H2S04  (1) 

An  excess  of  SO3  should  be  avoided  because  it  will  result  in  the  formation  of  (102)28207. 
The  sulfuric  acid  can  be  removed  by  washing  with  anhydrous  trifluoroacetic  acid  under 
anhydrous  conditions. 

During  the  course  of  this  study,  it  was  found  that  commercially  available  samples  of 
I2O5  were  not  suitable  for  the  neutron  experiments  because  they  contained  small  amounts 
of  hydrogen  containing  compounds.  These  not  further  identified  impurities  could  not  be 
removed  by  heating  the  samples  to  200  °C  in  a  vacuum  for  36  hours.  Pure  I2O5  suitable 
for  the  neutron  experiments  was  then  prepared  by  the  thermal  decomposition  of  I2O6: 

2I206  130tol70°c  >  2  I205  +  02  (2) 

This  is  an  exothermic  reaction  and  care  must  be  taken  to  remove  the  heat  of  reaction, 
especially  for  large  reaction  batches.  Otherwise,  the  reaction  mixture  can  overheat  which 
results  in  thermal  decomposition  of  I2O5  to  I2.  Such  an  event  was  observed  by  us  when  a 
250  mL  glass  flask  loaded  with  10  mmol  of  I20(,  was  heated  to  145  °C  in  an  oil  bath, 
resulting  in  a  run-away  exothermic  decomposition  reaction  resulting  in  partial 
decomposition  of  I2O5  to  I2. 

All  reactions  were  carried  out  under  dry  nitrogen  using  standard  Schlenk  techniques. 
Non-volatile  materials  were  handled  in  the  dry  nitrogen  atmosphere  of  a  glove  box. 
Glassware  was  heated  out  under  vacuum  before  use.  Raman  spectra  were  recorded 
directly  in  the  Teflon  reactors  in  the  range  4000-80  cm'1  on  a  Bruker  Equinox  55  FT-RA 
spectrophotometer,  using  a  Nd-YAG  laser  at  1064  nm  with  power  levels  of  200  mW. 
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Infrared  spectra  were  recorded  in  the  range  4000-400  cm"1  on  a  Midac,  M  Series,  FT-IR 
spectrometer  using  KBr  pellets.  The  pellets  were  prepared  inside  the  glove  box  using  an 
Econo  mini-press  (Barnes  Engineering  Co.)  and  transferred  inside  a  closed  container  to 
the  spectrometer  before  placing  them  quickly  into  the  sample  compartment,  which  was 
purged  with  dry  nitrogen  to  minimize  exposure  to  atmospheric  moisture  and  potential 
hydrolysis  of  the  sample.  The  starting  materials  HIO3,  H5IO6,  H2SO4  and  oleum  (all 
Aldrich)  were  used  without  further  purification.  Trifluoroacetic  acid  (SynQuest  Labs, 
Alachua,  FL)  was  freshly  destilled  from  P2O5  prior  to  use.  I2O6  was  prepared  by  a 
modified  literature  synthesis  [8,  9],  I2O5  was  prepared  by  thermal  decomposition  of  I2O6. 
The  purity  of  the  I2O6  and  I2O5  samples  was  checked  by  their  Raman  and  IR  spectra. 

Preparation  of  diiodine(V,VII)oxide,  I2O6 

Finely  ground  HIO3  (100.41  g,  0.571  mol)  and  H5IO6  (130.10  g,  0.571  mol)  were 
loaded  into  a  Schlenk  flask  and  cone,  sulfuric  acid  (400  mL)  was  added.  The  mixture  was 
heated  to  90  °C  and  stirred.  After  about  15  minutes,  a  yellow  suspension  has  formed.  The 
flask  was  now  removed  from  the  heat  source  and  30%  oleum  (184  mL)  was  added  under 
vigorous  stirring  to  the  hot  suspension.  The  mixture  was  stirred  at  ambient  temperature 
for  16  hours.  The  mixture  was  filtered  through  a  fine  porcelain  filter  frit  and  the  yellow 
solid  washed  ten  times  with  200  mL  of  dry  trifluoroacetic  acid.  The  solid  was  then  dried 
in  a  vacuum  at  ambient  temperature  for  14  hours.  We  thus  obtained  197.2  g  of  finely 
powdered,  yellow  I206  (weight  expected  for  0.571  mol  I2O6:  199.7  g). 

Preparation  of  diiodine(V)oxide,  I2O5 

I2O6  (2.50  g,  7.15  mmol)  was  loaded  into  a  glass  ampule  equipped  with  a  grease-less 
Kontes  HiVac  valve  with  a  Teflon  stopcock.  After  the  ampule  was  evacuated,  the  valve 
was  closed  and  the  vessel  placed  in  an  oil  bath  at  ambient  temperature.  The  bath  was  then 
heated  to  130  °C.  After  30  minutes,  the  temperature  was  raised  to  150  °C.  After  another 
30  minutes,  the  temperature  was  raised  to  160  °C  and  after  another  30  minutes,  raised  to 
170  °C.  The  solid  was  kept  at  this  temperature  for  12  hours  after  which  it  had  turned 
completely  white.  The  ampule  was  allowed  to  cool  to  ambient  temperature  and  connected 
to  a  glass  vacuum  line.  The  amount  of  non-condensible  gas  produced  inside  the  ampule 
was  determined  by  p,  V,  and  T  measurements  to  7.2  mmol).  We  thus  obtained  2.34  g  of 
finely  powdered,  white  I2O5  (weight  expected  for  7.15  mmol  I2O5:  2.39  g). 
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The  crystal  structures  of  anhydrous  I2O4  and  I2O5  over  the  temperature  range  of  4-300 
K  were  investigated  by  neutron  powder  diffraction  using  the  POWGEN  facility  at  the 
Spallation  Neutron  Source  (SNS)  of  Oak  Ridge  National  Laboratory.  Inelastic  neutron 
scattering  was  carried  out  at  low  temperature  (10  K)  using  the  chopper  spectrometer, 
ART,  also  at  SNS  [ref].  The  samples  were  sealed  inside  a  helium-filled  container 
throughout  the  experiments  to  avoid  interaction  with  the  ambient  atmosphere.  The 
extraordinary  large  neutron  incoherent  scattering  cross  section  from  hydrogen,  during  a 
long  run,  serves  as  a  sensitive  tester  of  the  presence  of  water  or  H-containing  impurities. 
We  find  no  evidence  of  incoherent-scattering  background  or  local-mode  vibrations  due  to 
hydrogen  species  in  the  samples. 

Owing  to  the  larger  size  and  more  polarizable  orbitals  of  the  I  atom  than  those  of  the 
O  atom,  the  binary  system  of  LOv  accommodates  the  increasing  oxygen  content 
structurally  via  adjustments  of  the  coordination  numbers  of  the  I  atoms  and  the  I-O-I 
bond  angles  but  always  stay  in  between  the  molecular-  and  framework-type  crystal 
structures.  Figures  1  (a)-(c)  illustrates  the  molecular  units  in  the  I2O4  (monoclinic, 
P2X  / c,  Z  =  4),  I2O5  (monoclinic,  P2l  / c,  Z  =  4),  and  I2O6  (triclinic,  PI,  Z  =  2)  in  the 
crystal.  In  all  cases,  neighboring  molecules  are  connected  by  some  intermolecular  1-0 
bonds  (2.05-2.9  A)  that  are  significantly  shorter  than  the  van  der  Waals  (vdW)  distance 

O 

(-3.5  A).  Furthermore,  each  structure  preferably  features  longer  intermolecular  1-0  bonds 

along  a  crystallographic  direction,  i.e.,  b ,  c ,  and  a  directions  for  I2O4,  I2O5,  and  I2O6, 
respectively.  H20  or  SO4  acid  molecules  tend  to  enter  the  lattice,  forming  inserted 
layered  structures  to  break  up  the  IOx  framework,  thereby  rendering  the  LOv  system 
thermodynamically  unstable. 

Computations  were  performed  using  the  Quantum  ESPRESSO  [10],  a  package  based 
on  the  density  functional  theory  (DFT),  plane  wave,  and  pseudopotential.  The 
pseudopotential  for  I  and  O  atoms  were  generated  by  the  method  of  Troullier  and  Martins 
[11]  in  conjunction  with  the  generalized-gradient  approximation  (GGA)  for  the 
exchange-correlation  functionals  [12].  To  account  for  the  VDW  interaction  between 
molecules,  an  empirical  VDW  correction  proposed  by  Grimme  was  incorporated  [13]. 
The  dispersion  correction  is  given  by 
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(3) 


■^disp  s6  5 /damply)’ 

. .  •  lia 

KJ  U 

where  Q  and  RtJ  denote  the  dispersion  coefficient  and  interatomic  distance  between  the 
ith  and  /h  atoms,  se  is  a  global  scaling  factor  that  depends  only  on  the  density  functional. 
A  damping  function,  /damp,  is  introduced  to  ensure  that  the  dispersion  correction  is 
negligible  for  small  Ry.  During  the  structural  optimization  using  the  scheme  of  damped 
variable-cell-shape  molecular  dynamics  (MD)  [14],  vdW  corrections  have  to  be 
introduced  to  bring  the  lattice  constants  matching  the  experimental  values.  The  final 
optimization  introduced  as  much  as  5-8%  corrections  for  the  b,  c,  and  a  lattice  constant  of 
I2O4,  I2O5  and  I2O6,  respectively,  and  about  one  third  of  these  values  for  other  lattice 
constants.  Specifically,  we  have  to  modify  the  global  scaling  factor  from  its  nominal 
value  of  S6  =  0.75,  to  1.0  and  2.0  respectively  for  I2O5  and  I2O6  in  order  to  account  for 
their  experimental  unit-cell  volumes.  This  anisotropy  results  from  the  tight  molecular 
netting  in  the  ac,  ab,  and  be  planes  in  the  I2O4, 12O5  and  I2O6  structures,  respectively,  see 
Figs.  1  (a)-(c).  This  characteristic  is  corroborated  by  experiments,  for  example,  through 
the  anisotropic  thermal  expansion  of  the  I2O6  lattice  favoring  the  a  direction  as  observed 
by  neutron  diffraction  (Fig.  1  (d)).  Due  to  the  large  difference  in  the  atomic  size  between 
I  and  O,  the  theoretical  implement  of  the  VDW  correction  has  a  non-trivial  impact  on  the 
bond  lengths  and  bond  angles.  In  general,  the  intermolecular  bond  lengths  agree  with 
experimental  values,  but  the  intramolecular  bond  lengths  are  overestimated  by  2-3%  and 
up  to  ~5%  in  some  cases. 

Dynamically,  the  effect  of  VDW  forces  correlates  not  only  with  the  variation  of  the 
vibrational  frequencies  but  also  more  rigorously  with  the  phonon  density  of  states 
(PDOS).  Therefore,  a  quantitative  comparison  of  the  combined  structural  and  dynamic 
data  with  first-principles  calculations  permits  a  systematic  assessment  of  the  interplay 
between  atomic/molecular  forces  and  the  structural/thermodynamic  properties.  The 
phonon  calculation  is  performed  by  a  diagonalization  of  the  dynamical  matrices  first 
computed  on  a  2x2x2  and  later  interpolated  to  a  6x6x6  mesh  in  the  first  Brillouin  zone 
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FIG.  1.  (a-c):  The  crystal  structures  of  ^Oy  (y  =  4,  5,  and  6),  where  the  unit  cells  are 
outlined  by  solid  lines  and  the  molecular  units,  I2O4,  I2O5,  and  I4O12,  are  identified  by 
atoms  denoted  in  solid  circles  (large  I,  small  O).  The  molecular  chain  along  the  re¬ 
direction  in  I2O4,  layer-like  substructure  in  the  ab- plane  in  I2O5,  and  framework 
configuration  in  I206  are  noticeable.  This  in  turn  favors  the  formation  of  longer  and 

weaker  intermolecular  1-0  bonds  (dashed  lines)  along  the  b ,  c ,  and  a  directions  for 
I2O4,  I2O5  and  I2O6,  respectively,  as  is  also  evident  as  anisotropic  vdW  interaction  in  the 
molecular-dynamics  simulation,  (d)  The  change  in  the  lattice  parameters  of  I2O6  relative 
to  those  at  300K,  showing  a  softer  c-direction. 


using  density  functional  perturbation  theory  [15]  and  the  visualization  of  the  lattice-, 
intermolecular-  and  intramolecular  vibrational  modes  by  the  Molden-4.7  software. 
Inelastic  neutron  scattering,  on  the  other  hand,  measures  the  neutron-weighted  (NW) 
PDOS  according  to 

2M/e2W(Q)  E  \  -  c-  • 

G(E)  =  — y\  ~ — 9 - -S(Q,E))DM  ^^(E)  (4) 

ft1  \  Q 2  n+ 1  /  . 
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where  c;,  a;,  M,-,  and  F£E)  are  the  concentration,  scattering  cross  section,  mass  and  partial 
PDOS,  respectively,  for  the  zth  atomic  species.  M  is  the  mean  sample  mass,  n  is  the  Bose 
population  factor,  and  S(Q,  E)  is  the  observed  neutron  scattering  function  for  one-phone 
excitations  as  a  function  of  energy  E  and  wavevector  Q.  The  Debye-Waller  factor 

e2lV(Q)  q2  fQr  one.pilonon  processes  can  be  estimated  from  diffraction  data.  (...}  and 

represents  the  average  over  a  wide  range  of  observed  Q  values.  The  experimental 
NWPDOS  is  compared  with  that  obtained  from  Eq.  (4)  using  the  F,(E)  from  MD 
simulations  convoluted  with  the  instrumental  resolution  function. 

The  calculated  and  measured  NWPDOS  for  I2O5  and  I2O6  as  well  as  the  experimental 
Raman  and  IR  frequencies  are  shown  in  Fig.  2.  Firstly,  the  lattice  modes  below  -15  meV 
extend  to  the  internal  molecular  vibrations  of  higher  energies  without  obvious  division, 
reflecting  the  comparable  bond  strengths  connecting  the  atoms  within  the  I2O4, 12O5,  and 
I4O12  units  and  the  atoms  between  neighbor  molecules.  For  I2O4  (not  shown)  and  I2O5, 
there  is  a  -10  meV  wide  gap  at  around  60  meV  dividing  the  bending/rocking  modes  and 
the  stretch  vibrations,  but  for  I2O6  such  gap  does  not  exist.  Moreover,  the  1 -phonon  cutoff 
is  seen  at  103,  109,  and  114  meV  for  I2O4,  I2O5  and  I2O6,  respectively.  These 
observations  are  consistent  with  the  stiffening  of  the  lattice  by  the  stronger  framework¬ 
like  connection  of  the  molecular  units  with  increasing  oxygen  content.  Secondly,  while 
the  Raman  and  IR  zone-center  modes  agree  well  with  the  neutron  peak  positions,  the 
calculated  NWPDOS  shows  some  marked  differences.  For  example,  For  I2O5  the 
calculation  predicted  higher  upper  frequencies  of  molecular  bending  modes  and  lower 
stretch  mode  energies  than  the  observed  values.  This  is  not  surprising  because  the  stretch 
modes  are  sensitive  to  the  bond  lengths;  their  low  values  are  the  direct  result  of  the 
overestimated  intramolecular  1-0  bond  lengths  resulted  from  the  vdW  correction. 
Fikewise,  the  vdW  correction  underestimated  the  short  intermolecular  1-0  bond  lengths 
hence  the  higher  predicted  frequencies  of  the  upper  bending  modes.  This  happens  also  in 
the  case  of  I2O4  through  an  examination  of  the  calculation  PDOS  and  the  Raman  and  IR 
frequencies  [11,  Ellestad  et  al.].  For  I206  the  calculated  1-0  band- stretching  frequencies 
are  clearly  underestimated  but  the  framework-like  crystal  structure  apparently  harden  the 
bending  modes  thereby  filling  the  gap  that  is  seen  in  I2O4  and  I2O5.  Finally,  the  observed 
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phonon  spectra  are  not  as  sharp  as  the  calculated  ones  perhaps  due  to  the  combination  of 
residual  multi-phonon  and  small-size  crystallite  broadening  effects. 


FIG  2.  The  measured  and  calculated  neutron- weighted  phonon  densities  of  states  for  I2O5 
and  I2O6  that  are  normalized  to  have  the  same  integrated  area,  and  the  Raman  and  IR 
data.  The  observed  IR  frequencies  about  50  meV  are  represented  by  tick  marks  and  the 
measured  Raman  spectra  are  displayed  in  arbitrary  units. 

The  MD  simulations  account  reasonably  well  the  crystal  structures  and  lattice 
dynamics  of  the  kCf  system  notwithstanding  the  10-15%  energy  shifts  of  the  high-energy 
modes  caused  by  the  compromise  between  inter-  and  intramolecular  interactions  from  the 
VDW  correction.  We  have  calculated  the  equation  of  state,  thermal  expansion,  heat 
capacity,  adiabatic  bulk  modulus  and  thermal  Grtineisen  parameters  (the  weighted 
average  of  the  mode  Griineisen  parameters)  for  I2O4,  I2O5  and  I2O6.  Figure  3  shows  the 
effects  on  the  calculated  equations  of  state:  the  substantial  overestimation  without  the 
vdW  correction  and  the  adjustment  of  the  s6  parameter,  and  the  improvement  by 
including  the  thermal  contribution  of  phonons  under  the  quasi-harmonic  approximation. 
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The  contribution  of  thermal  vibrations  (including  the  zero-point  motion)  to  the 
equilibrium  volume  as  a  function  of  pressure  is  evident  and  amounts  to  about  0.9  to  2.2%. 
Interestingly,  these  values  are  comparable  to  those  of  MgO  and  Mg2Si04  [16]  in  spite  of 
the  large  difference  in  the  bulk  modulus  between  the  latter  3D  framework  crystals  (-170 
GPa)  and  the  iodine  oxides  (-30  GPa).  The  three  iodine  oxides  have  similar  thermal 
pressure  gradient  of  -1.7  GP/1000K.  As  a  result,  the  magnitude  of  the  bulk  modulus, 
i,o  >  i,o  >  w,,foU°ws  the  order  ofthc density,  >  ,0  >  ,0. 


FIG  3.  The  calculated  equations  of  state  for  I2O4  (lines),  I2O5  (lines  affixed  with 
diamonds),  and  LOf,  (lines  affixed  with  circles).  For  each  I2Oy  member  calculations 
without  vdW  force  corrections  (dashed  lines),  with  vdW  force  correction  but  not 
including  thermal  contribution  from  phonons  (dotted  lines),  and  with  both  vdW  force 
correction  and  phonon  contributions  at  300K  (solid  lines)  are  shown.  Experimental  cell 
volumes  for  4  formula  units  at  ambient  pressure  are  shown  by  large  diamonds. 

In  summary,  neutron,  Raman  and  IR  measurements  of  the  crystal  structures, 
vibrational  mode  frequencies,  and  phonon  densities  of  states  of  newly  synthesized  pure, 
anhydrous  I2O5  and  LOft  samples  provided  the  crucial  experimental  data  for  the 
optimization  of  both  the  structures  and  dynamics  of  I2O4, 12O5,  and  LO^  via  ab  initio  MD 
simulations.  Important  thermodynamic  properties  were  computed  for  further 
understanding  of  the  materials  under  extreme  conditions. 
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